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Abstract 

Survival of floret primordia initiated seems critical for tlie determination of grain number and yield in wheat, and 
understanding what determines floret mortality would help in the development of more robust physiological models 
of yield determination. The growth of the juvenile spiltes has been frequently considered the determinant of grain 
number, implying that floret development would depend on resource availability and that the onset of floret death 
would be related to spike growth. However, this model has been recently challenged from a study concluding that 
floret death started when the most advanced floret primordia reached a particular developmental stage. As the few 
previous studies on this relationship involved photoperiod treatments which affect both floret development and the 
onset of spike growth, conclusions cannot be considered mechanistic. This comprehensive study analysed in detail 
floret development in wheat as affected by resource availability (mainly soil nitrogen levels) and found that the onset 
of floret death may occur when development of the most advanced florets ranged from stages 5 to 9 and that the 
average and standard deviation of floret developmental stage coinciding with the onset of floret death was not related 
to the level of availability of resources. These results provide further support to the model relating the onset of floret 
death with the initiation of active growth of the juvenile spike in which florets are developing. 
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Introduction 

To further increase wheat yield, it may be useful to under- 
stand the mechanisms controning grain number determi- 
nation (Slafer, 2003; Fischer, 2011; Reynolds et a!., 2012). 
Grain number is related to the number of fertile florets 
(Kirby, 1988; Slafer and Andrade, 1993; Miralles et al., 
1998, 2000) and the latter is the end result of the dynam- 
ics of floret primordia generation/degeneration during the 
stem elongation phase (Kirby, 1988). Most frequently, what 
determines the number of fertile florets is the survival of pri- 
mordia initiated rather than the number of primordia pro- 
duced (Ferrante et al., 2013 and references quoted therein). 
Understanding what determines floret mortahty might be 
helpful to have a better physiological model of the process 
(Slafer et a/., 2005). 



As it has been shown several times that the numbers of 
both fertile florets and grains are related to spike dry weight 
at anthesis (Miralles et al., 1998; Acreche et al., 2008; Fischer, 
201 1; Gonzalez et al., 201 1; Ferrante et al., 2012), it would be 
reasonable to assume that floret development would depend on 
resource availability. In fact, there seemed to be a relationship 
between the fate of individual florets, determining floret sur- 
vival, and spike growth dynamics as affected by photoperiod 
manipulation during stem elongation (Gonzalez et al., 2003b, 
2005; Serrago et al., 2008). The introgression of Rht genes did 
also seem to have increased grain number through increasing 
floret survival related to higher rates of spike growth before 
anthesis (Miralles et al., 1998). This model, assuming that the 
number of fertile florets is mainly the outcome of a simple 
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trophic process, would be further supported if it is determined 
that what triggers floret mortahty is also related to spike 
growth. In principle, and in line with the model of floret sur- 
vival, it has been proposed that the onset of floret death would 
not take place before the beginning of active spike growth 
(Fig. 1, left). A tough problem to directly test the assumption 
is that determining the dynamics of floret initiation and death 
to assess the onset of floret death requires a huge amount 
of work (dissecting very often individual floret primordia to 
determine their developmental stage) and therefore most com- 
monly the few results which are available are derived from the 
imposition of only few treatments (and the lower the number 
of treatments, the higher the uncertainty on the strength of the 
conclusions based on coordination, or lack of, between pro- 
cesses). The proposition that the onset of active spike growth 
would trigger floret death was based in a couple of studies by 
Gonzalez et al. (2003b, 2005), in which plants were subjected 
to contrasting photoperiods during stem elongation. On the 
other hand, it has been recently concluded that the onset of 
floret death would be a 'pure' developmental process, unre- 
lated to spike growth, governed by the stage of development 
of the most advanced floret primordium of the spike, the most 
proximal floret (Fl) of central spikelets (Bancal, 2008, 2009; 
Fig. 1, right). Again, this conclusion was based on a rather 
small database (only three environmental treatments: a con- 
trol, an extended photoperiod, and shading of 28% imposed 
on six cultivars). More recently, Gonzalez et al. (2011) reana- 
lysed a substantially larger database, in which timing of spike 
growth and development of florets were mainly, though not 
only, affected by photoperiod treatments, and concluded that 
onset of floret death was most likely related to the onset of 
spike growth at maximum rate. However, due to the nature of 
the treatments imposed, there was also a correlation between 
the onset of spike growth and floret developmental rates. As 
longer photoperiods advance both spike growth and floret 
development, resolving the issue would require treatments 
affecting growth but not development. 

This study undertook a number of experiments in which 
(i) wheat plants were exposed to treatments affecting crop 
growth but not crop development (mainly nitrogen fertiliza- 
tions) and (ii) the developmental stage of floret primordia 
determined in detail at different positions in the spikes during 
stem elongation. This study aimed to determine whether the 
onset of floret death was related to the developmental stage 
of Fl in central spikelets by reanalysing the large database 
based on these experiments. 



Materials and methods 

General conditions 

Six experiments were conducted, during four consecutive grow- 
ing seasons from 2006-07 to 2009-10 (Table 1). Four experiments 
(experiments 1, 2, 3, and 6) were carried out outdoors in micro- 
crops, within large cubic containers ( 1 m lieiglit and 1 x 1 m^ surface; 
details available in Ferrante et a I., 2010) at the School of Agronomy, 
University of Lleida, Spain (41° 37' 50" N 0° 35" E; altitude 180 m) 
while the other two were field experiments sown at Gimenells (41° 
39' 11" N 0° 23' 28" E; altitude 258 m). In order to ensure a low 



nitrogen (N) availability in experiments 1,2, 3, and 6, the contain- 
ers were filled with a sand/soil mix (3:1 v/v). In each experimental 
unit, mineral N availability at sowing was equivalent to 70 kg N 
ha"' (experiment 1), 30 kg N ha"' (experiment 2), and 20 kg N lia"' 
(experiments 3 and 6). Before sowing, phosphorus fertilizer (tri- 
ple superphosphate, 20kg P ha"') was uniformly mixed within the 
upper 20 cm of the soil mixture. Before tillering, micronutrients were 
added. Experiments 4 and 5 were conducted under field conditions 
in a soil classified as calcixerept petrocalcic (FAO, 1990). The main 
environmental differences between the two experiments was the soil- 
N contents and irrigation regime, although they also differed more 
slightly in sowing dates and soil depth (0.75 m in experiment 4 and 
more than 1 m in experiment 5). 

All experiments were sown within the normal sowing dates of the 
region and at normal sowing densities (Table 1 ). Microcrops (experi- 
ments 1, 2, 3, and 6) were sown with a dedicated system to warrant 
perfectly uniformly distributed seedlings (as described and illus- 
trated in Ferrante et al., 2010). 

In all experiments, diseases and insects were prevented or con- 
trolled by spraying recommended fungicides and insecticides at the 
doses suggested by their manufacturers. In addition, weeds were 
removed by hand (microcrops) or controlled by spraying selective 
herbicides (field). 

The mean daily global radiation (M J m"^), mean temperature (°C), 
and accumulated precipitation (mm) were recorded daily in meteor- 
ological stations belonging to the Catalonian Agro-meteorological 
network, located close to the experimental sites, and are shown for 
the major developmental phases in each experiment in Table 2. 



Treatments and experimental designs 

Although treatments varied among experiments, in all there were 
different levels of N availabilities and/or different cultivars, in some 
cases combined with contrasting water regimes (Table 1). 

During 2006-07 and 2007-08, cultivar Claudio (modern and well- 
adapted durum wheat) was grown in micro-crops with contrasting 
N (experiments 1 and 2) and water (experiment 2) availabilities 
(Table 1 ). In the third growing season (2008-09), the main treatments 
in the three experiments were four modern (released between 1987 
and 2002) semidwarf durum wheat cultivars: Claudio, Donduro, 
Simeto, and Vitron (Table 1). In the experiment in large containers 
they were subjected to different N availabilities (experiment 3). The 
other two were field experiments with contrasting conditions: rain- 
fed (experiment 4) and irrigated (experiment 5). In the fourth grow- 
ing season (2009-10; experiment 6), two of these cultivars, Donduro 
and Vitron (selected because they were the two highest yielding cul- 
tivars in the third year, while exhibiting contrasting response to N 
in terms of spike fertihty; Ferrante et al., 2012), were exposed to a 
factorial combination of two contrasting N and water availabilities. 

Judged by the amount of certified seed produced, Claudio was the 
most widely grown cultivar in 2009/10 and is particularly adapted 
to the region; therefore, it used as a check for testing new cultivars 
in the northern growing area of Spain. Vitron and Simeto are cul- 
tivars of much wider adaptation than Claudio as they are checks 
not only for the northern but also the southern region of Spain 
(GENVCE, Group for the Evaluation of the New Cereal Varieties 
in Spain; Anonymous, 2008). Donduro is the newest cultivar and 
soon after its release it showed in comparative yield trials to perform 
consistently better than Vitron and Simeto over the wide range of 
conditions explored in that network (Anonymous, 2008). The four 
cultivars chosen are all modern and well adapted to the growing 
conditions and therefore their general phenological patterns are very 
similar and their plant heights are within the optimum maximizing 
yield (Ferrante et al., 2012). 

Measurements and analyses 

From terminal spikelet initiation to flowering, one plant per experi- 
mental unit was selected at random and harvested twice or thrice a 



Onset of floret death and wheat floret development | 2861 
Alternative models for the onset of floret death 



Trophic model proposed by 
Gonzalez et al. (2005; 2011) 



Pure developmental model 

proposed by Bancal (2008; 2009) 
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Fig. 1. Diagram schematically describing two alternative models for processes and signals responsible for triggering the onset of floret 
death (downward arrows in the top panel): the trophic model (left side) in which floret death would be triggered by the dynamics of 
spike dry weight (SDW) between terminal spikelet initiation (TSI) and anthesis (A), and the 'pure' developmental model (right side) in 
which floret death would be triggered by the developmental stage of the most advanced floret primordium. Thus, floret death starts 
either at the onset of rapid spike growth in the trophic model, or when F1 of central spikelets reach the stage 8 of the Waddington ef a/. 
(1 983) scale in the developmental model. Support for the trophic model comes from studies in which by reducing the length of the 
TSI-A period the onsets of both spike growth was floret death were advanced in parallel. Support for the developmental model comes 
from studies in which the onset of floret death was correlated with the timing required for F1 to achieve the stage 8 in the mentioned 
scale. The bottom part of the scheme represents the model of floret survival from onset of floret death to grain number, where there is 
consensus in the literature on the paramount relevance of floret survival in determining grain number in wheat. 
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week. The spike of the main shoot was dissected under a binocu- 
lar microscope (Leica MZ 7.5, Leica Microsystems, Heerbrugg, 
Switzerland) to count the total number of floret primordia. In addi- 
tion, at each sample the developmental stage (floral score) of each 
floret within particular spikelets was determined following the scale 
of Waddington et al. (1983), as described in Ferrante et al. (2010, 
2013). The spikelets analysed were those in basal (fourth spikelet 
from the base), central (middle spikelet), and apical (fourth spike- 
let from the top) positions of the spike. Florets within the spikelets 
were numbered from Fl (floret primordium closest to the rachis) to 
Fn (the most distal floret primordium of the particular spikelet, this 
being normally F8). Floret primordia were considered as fertile flo- 
rets when they were at W 10 or immediately before that stage (when 
the stigmatic branches were curved with green anthers). 

To determine the effect of treatments on the different variables, 
the data were subjected to analysis of variance and the relation- 
ships between variables were determined by regression analysis 
(SAS, 2002). Developmental score of Fl at the onset of floret death 
was determined, for each combination of resources availability and 
growing seasons. For this purpose, the number of living floret pri- 
mordia was plotted against the stage of development of Fl, and the 
stage of Fl at the onset of floret death was determined. For ana- 
lysing all data from different spikelets and treatments together, the 
number of living floret primordia were standardized as a proportion 
of the maximum number of floret primordia initiated in each par- 
ticular cases (Gonzalez et al., 201 1). 

During the last growing season (2009-10), spike length and spike 
dry weight were also determined in each of the samples to describe 
floret development. The dynamics of spike length and dry weight 
were determined by fitting the data with a bilinear model and the 
onset of rapid spike growth was determined. 
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As expected, N did not affect the developmental rate of flo- 
rets which are normally fertile at anthesis. This is illustrated 
in detail for the most proximal floret (Fl) in central spike- 
lets under low and high N availability (Fig. 2, upper pan- 
els). Development of Fl was never affected significantly by 
N availability (Supplementary Table SI, available at JXB 
online). Development of F2, which is also most frequently 
a fertile floret even under low N levels, did not respond to N 
either (data not shown). For more distal florets, like F4, which 
under some conditions becomes a fertile floret but under 
other conditions does not survive and dies, the effect of N on 
their development can be seen in some cases (Supplementary 
Table SI), mainly those in which F4 was fertile under high 
N but died at earlier stages under low N availability, as illus- 
trated for F4 in experiment 1 (Fig. 2, bottom panels). 

As N did not affect the developmental pattern of the most 
proximal floret in central spikelets, it proved to be the type of 
treatment required to test whether the onset of floret death is 
actually determined by the timing of achieving a certain stage 
of development in these florets. If this is so, N availability (and 
the other treatments evaluated in the present study. Table 1) 
should have not affected the timing of the onset of floret death. 
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The dynamics of the relative number of living florets through 
time was very similar between all experimental years and 
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Fig. 2. Relationship between floret developmental score (scale of Waddington et al., 1 983) and thermal time from sowing of either 
the most proximal (and most advanced) floret primordia (F1) or the fourth floret primordia from the rachis, a relatively distal floret (F4) in 
central spikelets under contrasting N treatments in experiment 1 . 



spikelet positions. In all cases, the number of floret primor- 
dia increased sharply, reached a maximum, and then declined 
markedly to end up in a particular number of fertile florets, 
which normally ranged between 20 and 40% of the maximum 
number of floret primordia initiated. This process of floret 
generation/degeneration occurred while the proximal floret of 
the central spikelet progressed in its development (Fig. 3, top 
panel). 

Regardless of the noticeable similar general pattern, there 
was not a clear single developmental stage of the most prox- 
imal floret of central spikelet synchronic with the onset of 
floret death: depending on the particular case (experiment, 
treatments, spikelet positions), the onset of floret death coin- 
cided with a developmental score of the most proximal flo- 
ret of central spikelets, ranging from 5 to 9 (Fig. 3, bottom 
panel). 

In order to see whether the variation was induced by 
the different environmental conditions explored, the data 
were reanalysed within particular environments. For all, 
a strong association was found between the onset of floret 



death and the stage of development of the most advanced 
floret primordia in the spikes (Fig. 4, top panel). Again, 
within environmental conditions, floret death started when 
the developmental score of Fl of central spikelets ranged 
roughly from 5 to 9, depending on the cultivar and spike- 
let position considered (Fig. 4, bottom panels). In addition, 
similar to the finding using the developmental score of Fl in 
central spikelets, there was not a clear single floret develop- 
ment score associated with the onset of floret death when 
considered the most proximal floret of the corresponding 
spikelets (basal, central, and apical spikelet positions), either 
analysing all data pooled together (Supplementary Fig. SI) 
or analysing data within particular environmental condi- 
tions (Supplementary Fig. S2). 

The developmental score of Fl of central spikelets 
that resulted most frequently related to the onset of flo- 
ret death varied across experimental conditions, between 
6.5 under low resources at the University of Lleida and 
8 under high inputs in Gimenells (Fig. 4, bottom panels). 
However, the variation in stage of development of Fl at 
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Fig. 3. (Top) Relationship between tine number of living florets, 
in relative terms to standardize for differences in maximum 
number of floret primordia depending on treatments and 
spikelet positions, and score of floret F1 in central spikelets. 
Open circles corresponding to the dynamics of all data (in each 
particular case unidentified in this figure always the values peak 
to 1 before the onset of floret death). To illustrate the range 
of differences in score of F1 in central spikelets at the onset 
of floret death in the wide range of conditions analysed, we 
identified the two extreme cases in which the onset of floret 
death coincided with F1 in central spikelets being at stages as 
early as a score of 5 (closed triangles) or as late as a score of 9 
(closed rhomb). (Bottom) Frequency of developmental scores of 



which floret death started did not seem to be consistently 
affected by the particular growing conditions (Fig. 5a and 
Supplementary Fig. S3a). To analyse this more quantita- 
tively, assuming the yield broadly reflects the availability 
of resources, the mean score of Fl in central spikelets at 
the onset of floret death was related with the average yield 
in each of the conditions: there was no consistent trend 
(Fig. 5b and Supplementary Fig. S3b). There was no trend 
between the standard deviation in score of Fl at onset of 
floret death and yield in the particular conditions (Fig. 5c 
and Supplementary Fig. S3c), implying that the range in 
developmental score of Fl at onset of floret death was 
similar for the dramatically different growing conditions 
explored in this study. 

Spike length and spike dry weight 

Although spike length and spike dry weight are acceptably 
related when both are determined at flowering, the dynam- 
ics of spike length and dry matter accumulation were 
rather different (Fig. 6). The onset of rapid spike lengthen- 
ing took place at about 200 °C d before the onset of rapid 
spike dry matter accumulation (Fig. 6, top and middle 
panels). Thus, spike length can hardly be used as a proxy 
of spike dry matter if the aim is to relate early spike weight 
with dynamics of floret development. The only case in 
which spike dry weight was determined frequently between 
terminal spikelet and anthesis showed that the onset of flo- 
ret death occurred soon after the onset of rapid dry weight 
accumulation on the juvenile spike (Fig. 6). In fact, both in 
Donduro and Vitron for the last experiment, the onset of 
floret death followed the onset of spike growth though due 
to their difference in development it was coincident with 
different stages of Fl development in central spikelets (W 
8.2 and W 8.8 in Donduro and Vitron, respectively; Fig. 6, 
bottom panel). 

Discussion 

The main goal of this paper was to determine whether the 
onset of floret death was triggered by the achievement of a 
certain developmental stage (c.8 in the scale of Waddington 
et al., 1983), by the most advanced floret primordia (Fl of 
central spikelets) and independently of the growth of the 
spike. This was a conclusion derived from a study reported 
by Bancal (2009), implying that the onset of floret death 
would be, therefore, independent of spike growth dur- 
ing pre-anthesis (reinforcing a previous paper in which he 
inferred that spike growth may not determine number of 
fertile florets or grains; Bancal, 2008). If this were so, much 



the proximal florets (Fl) of central spikelets at which floret death 
started considering the whole dataset analysed. Bars represent 
standard errors. Pictures show the appearance of the primordia 
of F1 in central spikelets corresponding to two extreme and one 
intermediate conditions. 
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Fig. 4. (Top) Number of living florets (relative to the maximum in each case) plotted against the developmental stage of the most 
advance floret (F1) in central spikelets for each of the three spikelet positions considered (upper panels) analysed for each of the six 
environmental conditions explored throughout the studies analysed independently. (Bottom) Frequency of developmental scores of the 
F1 of central spikelets at which the onset of floret death took place in each environmental conditions studied. 



of the most accepted physiological models for the determi- 2011), would need to be conceptually revised, and ideas 
nation of grain number, based on interplay between spike on how manipulating genetically spike growth to maxi- 
growth and floret survival (e.g. Slafer et al., 2005; Fischer, mize grain number in wheat (Miralles and Slafer, 2007; 
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Fig. 5. (a) Boxplot of the stages of development of the most advanced floret (F1 ) in central spikelets at the onset of floret death for 
different experimental conditions (experiments and treatments). The horizontal line within each box is the median, and the bottom and 
top of each box represent the 25th and 75th percentiles, respectively. Crosses correspond to the mean. The whiskers represent the 
10th and 90th percentiles, respectively and the values outside this range are individual outliers (closed circles), (b and c) Relationships 
between the average (b) and the standard deviation (c) of the scores of floret F1 in central spikelets at the onset of floret death in each 
environmental condition analysed and yield achieved by the crop in those conditions (as indicator of resource availability). Environmental 
conditions were (i) high- and low-nitrogen levels in the microcrops (closed and open symbols, respectively) under irrigated (circles) or 
rainfed conditions (triangles), and (ii) high- or low-inputs (closed and open squares, respectively) in the field experiments. 



Foulkes et al., 2011; Reynolds et al, 2012) might be more 
doubtful. Therefore, testing the robustness of the relation- 
ship between the onset of floret death and achievement of 
a certain developmental stage in Fl of central spikelets 
is rather relevant. The mentioned studies by Bancal are 
extremely useful because of the detail provided for floret 
development under contrasting conditions but there are 
two major inconveniences to allow accepting their conclu- 
sions straightforwardly. First, the study included only very 
few environmental conditions to yield easily extrapolatable 
conclusions. Second, the likelihood of the onset of floret 
death be triggered by spike growth was discarded assum- 
ing that floret death should have started later than the tim- 
ing when the rate of spike growth is maximized, which in 



the analysis occurred when spike dry weight was already c. 
0.5 of the weight at anthesis (clearly much later than the 
actual onset of spike growth; see also discussion on this in 
Gonzalez et al., 201 1). 

The present study comprehensively tested the existence of 
such developmental triggering of floret death with a com- 
prehensive database of both dynamics of floret generation/ 
degeneration (from where the onset of floret death could be 
determined) and detailed developmental progress of partic- 
ular floret positions from terminal spikelet initiation to few 
days after anthesis, as affected by a number of treatments 
changing resource availability (mainly, but not only, soil- 
N availability levels), which are well recognized to greatly 
affect crop growth but not development. In agreement with 
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Fig. 6. Dynamics of the number of living floret primordia per spiketet 
(NLF, closed squares), spike lengtln (SL, inverted grey triangles), 
and spike dry weiglit (SDW, open triangles) from terminal spikelet 
initiation to antlnesis for Donduro (top panel) and Vitron (middle 
panel) under the irrigated and fertilized treatment of the last growing 
season. The bottom panel shows the relationship between the 
number of living floret primordia per spiketet and the developmental 
score of F1 in central spikelets for the same treatments. Dotted lines 
indicate the stage of F1 development at the onset of floret death. 



Gonzalez et al. (2011), who also analysed this relationship 
(but with a database mainly produced by manipulations of 
photoperiod), the current study found the onset of floret 
death to be rather independent of the developmental stage of 
the most advanced floret. In the database, floret death might 
start when Fl of central spikelets are anywhere between the 
stages 5 and 9 of the scale of Waddington et al. (1983); which 
might represent a difference in time of more than 200 °C 
d. The results of this study do therefore conflict severely with 
the conclusion that the onset of floret death may be simply 
a developmental process, giving further support then to the 
more-established model that the processes related to floret 
survival to determine number of fertile florets and grains are 
tightly related to the growth of the juvenile spikes during the 
stem elongation phase. 



Supplementary material 

Supplementary data are available at JXB online. 

Supplementary Fig. SI. Relationships and frequencies 
equivalents to those shown in Fig. 3 but corresponding to 
bottom, central, and apical spikelets. 

Supplementary Fig. S2. Number of living florets and fre- 
quencies equivalents to those shown in Fig. 4 but correspond- 
ing to bottom, central, and apical spikelets. 

Supplementary Fig. S3. Boxplot and relationships equiva- 
lents to those shown in Fig. 5 but corresponding to bottom, 
central, and apical spikelets. 

Supplementary Table SI. Spike dry weight at anthesis and 
parameters of the bilinear regression analysis of floret devel- 
opmental stage vs. thermal time for the most proximal and 
one distal florets of central spikelets of wheat grown under a 
wide range of conditions. 
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